Spectral Effects of Fast Response Cross Kerr Non-Linearity on Quantum Gate 
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According to idealized models, a strong Kerr non-linearity may be used to build optical quantum 
gates for optical quantum information processing by inducing conditional phase shifts on quantum 
states. Recently, Shapiro (PRA 73, 062305 (2006)) argued that for a Kerr medium with non- 
instantaneous but fast response, essentially no phase shift is induced on two-single-photon input 
states, and thus a quantum gate built from such a medium cannot work. Here we show that a fast 
response Kerr medium induces some but very little phase shifts on a two-single-photon input state, 
and it is insufficient for high fidelity quantum computation. We point out that this is caused by 
the medium imparting spectral entanglement to the input photons. We further show that a way to 
circumvent this problem and achieve a high fidelity gate, is to engineer the dispersion properties of 
the medium to give a dominant spectral effect over the non-instantaneous response, in addition to 
satisfying a phase matching condition. 

PACS numbers: 03.67.Lx, 42.50.-p 
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I. INTRODUCTION 

Research in quantum information sciences has shown 
that a large scale quantum computer has the ability 
to factorize large numbers, search through unsortcd 
databases and perform quantum simulations, more ef- 
ficiently than classical computers. In order to achieve a 
universal set of quantum logic operations, quantum en- 
tangling gates between qubits are required. Single pho- 
tons have been demonstrated to be good carriers of quan- 
tum information. However, photons interacts with each 
other weakly and this provides a major obstacle in build- 
ing an optical quantum computer. Nevertheless, various 
schemes have been proposed to construct a two-qubit 
quantum entangling gate. Linear optics with measure- 
ment induced non-linearity can be used to build an entan- 
gling gate but it suffers from a low success rate, which in 
other words means lots of resource are required to achieve 
high efficiency P, 0, S| ■ This has led to renewed interests 
in building a quantum gate with non-linear optics, which 
can provide stronger interaction between photons [4] _5| . 
Recent developments in non-linear optical quantum gates 
include gates such as the quantum Zeno gate [|| 0| and 
the Kerr non-linearity gate [1, Q . 



Cross-Kerr Quantum Gate 



can induce a it phase shift on the two-single-photon state 
|l)p|l) s without altering the |0) p |0) s vacuum state, and 
the |0)p|l) s and |l) p |0) s single-photon states. For |0)j 
being logical zero and \ l)j being logical one, this is equiv- 
alent to a single rail logic controllcd-sign gate. We note 
that since each mode may have at most one photon in 
our calculation, we do not include a self-phase modula- 
tion term in the Hamiltonian, because it will have zero 
contribution. For this Hamiltonian, the unitary operator 

is simply U = cxp • A n phase shift on the |11) state 

can be obtained if ^ = 7r, and this gives the aforemen- 
tioned controllcd-sign operation. If we encode our qubits 
in polarization encoding, say \H) is logical and \V) is 
logical 1, then this operation can be implemented as a 
gate in dual-rail logic , as shown in figure [TJ Along 
with arbitrary single-qubit operations that can be done 
with wavcplates , universal quantum computation can 
be achieved. 
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The Hamiltonian for a single frequency cross-Kerr 
medium with modes p and s can be written as: 

H = xa\a p a\a s (1) 

where \ denotes the interaction strength. With a suf- 
ficiently strong cross-Kerr interaction, the Kerr medium 



FIG. 1: Schematic of the cross-Kerr controlled-sign gate in 
dual-rail logic. PBS stands for polarization beam splitters. 



Since the response of a non-linear medium is frequency 
dependent, it is therefore important to know how its 
spectral response affects the spectral profile of photonic 
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input states, as well as the fidelity of the gate. Recently, 
Shapiro fl2j has argued that for a Kerr medium with 
non-instantaneous but fast response, there is essentially 
no phase shift induced on the two-single-photon input 
states. Thus quantum gates based on such media cannot 
work. Furthermore, for a Kerr medium with slow 
response, significant phase shift can be imparted to the 
state but the presence of phase noise in this regime pre- 
vents a high fidelity operation. Here we show that a fast 
response Kerr medium does induce some phase shifts on 
two-single-photon input states but insufficient for high 
fidelity quantum computation. We identify the problem 
as predominately arising from spectral entanglement, 
that is, we show that for separable two-single-photon 
input states, the output state from such a medium is 
inevitably spectrally entangled. We then further show 
that a way to circumvent the problem is to engineer the 
dispersion properties of the medium to give a dominant 
spectral effect over the non-instantaneous response, in 
addition to satisfying a phase matching condition. 

This paper is arranged in the following way. The next 
section discusses how we model the spectral properties of 
a cross-Kerr non-linear medium that has a fast but 
non-instantaneous response. We then derive the output 
state for a separable two photon input state and find 
the fidelity of the gate. In Section III, we indicate how 
the problem of spectral entanglement can be eliminated 
by introducing dispersion into the model and deduce a 
phase matching condition that can enhance the fidelity. 
We conclude in section IV. 



II. MODELLING THE SPECTRAL RESPONSE 
OF A KERR MEDIUM WITH FAST RESPONSE 



H(t) =x i3) e J v dr 3 (E p \r,t)E p (r,t)(^ dr ff (r) 
E s \r,t - r)E s (r,t- r) + Km p (r,r) 
+ ( J dT 9 (T')E p \v, t - r')E p (r, t - r') 

+ K m s {T / 7 r))E s \r 7 t)E s (r 7 t)) (3) 



Hereafter, we shall limit the analysis to one spatial 
dimension, the propagation direction. The expression for 
the electric field operator of mode j is: 



where Aj (u)j 



(4) 



hujj 



and S is the cross section 



area of the beam and rij{u>j) is the refractive index for 
mode j. We assume that Aj{wj) = Aj is slowly varying 
for the frequencies of interest, allowing it to be factored 
outside the integral. In this section, we are interested 
in the effects of non-instantaneous response, and so 
we assume that it dominates over the effects of dis- 
persion and hence we assume wavenumbcr kj is constant. 

The m(r, z) terms in the Hamiltonian are noise opera- 
tors that are required to preserve the commutation of the 
field operators when the medium has a non-instantaneous 
response 12j, [l5|, Il6|. In the case of self-phase modula- 
tion, Boivin ct al [16| showed that the noise term can 
be modelled as a collection of localized and independent 
harmonic oscillators with the expression: 



The process of Kerr non-linearity can be studied in the 
interaction picture, with the unitary evolution of a state 
vector given by: 



Ufa, i )|V) =exp fr{ - 
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m}) IV') (2) 



where T is the time ordering operator. Since the inter- 
action Hamiltonian for a cross-Kerr non-linearity does 
not necessarily commute at different times, therefore the 
Dyson series 13 , 3] shall be used to calculate the higher 
order terms of the unitary expansion. In the case of a 
time-commutable interaction Hamiltonian, the time or- 
dering operator has no effect and may be dropped, re- 
sulting in the usual Taylor series for the unitary expan- 
sion. The interaction Hamiltonian for a cross-Kerr non- 
linear medium with non-instantaneous response and no 
self Kerr effect, has the expression as shown in equa- 
tion ©. 



m(T,z)= I dLud^-dl{z)e lWT + H.C. (5) 



where G{lo) = J dTsin(uiT)g(T). Furthermore, in the case 
of under-damping, i.e. < T < 2Q, the response function 
for t > is 



9(r) 



KQ 2 



e.- rT/2 sin (V^ 2 - r 2 /4r) (6) 



where T is the damping coefficient and O is the resonance 
frequency of the medium. Shapiro [l2| commented that 
this model of noise also applies to cross-phase modula- 
tion. If the medium has a very high resonance frequency 
and a large damping coefficient, g(r) resembles a delta 
function, which gives the fast response limit that we 
are interested in. To illustrate this, for instance, lets 
consider the limit of critical damping, T = 251, and 
such that g{r) = KQ 2 Te~ nT . The maximum of this 
function is Kfle^ 1 and occurs at r = Furthermore, 
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the full width half maximum is proportional to 1/fi. 
Hence, for large f2, the response function resembles a 
delta function. For mathematical convenience in our 
calculation, in the fast response regime, we approximate 
g(r) « exp(— M 2 r 2 ) as a narrow gaussian function 
in time with r > and large Af . Here, M oc f2 and 
that in the limit where M — > oo, we have g(j) — » <5(t). 
As the response function tends to a delta function, 
equation §5§ shows that the noise term tends to zero 
and therefore in the fast response regime, the noise 
term has negligible contribution and the noise term 
in the Hamiltonian shown in equation (|3|) can be ignored. 

For the Hamiltonian, after integrating over r by as- 
suming large M, and similarly for the s mode, in the 
expression we have the term: 



enhance the gate fidelity Hence, for the purpose of illus- 
trating the spectral entanglement problem, it is sufficient 
to use the Taylor series, even though formally we should 
use the Dyson scries. Nonetheless, in the appendix, 
we will further discuss the validity of this approximation. 

Integrating the Hamiltonian in equation ([8]) over time 
t with bounds taken from — oo to oo by considering far 
field limits (HI, results in a factor of S(Auj). When we 
apply the energy conservation condition by integrating 
d(Auj), the Hamiltonian becomes 



H = H(t)dt 

J — OO 

=X L j j J dLd p+ dLd p -duj s+ al(uj p+ )a p (uj p -)al(uj s+ ) 
d s (uj p+ - + uj s+ ) cxp (~ ^ P+ 4Af 2 P ~^ ) ( 9 ) 



Thus the interaction Hamiltonian has the expression: 
XL 



H(t) 



dcu p + dtUp- duj s -\-doj s - 



' — oo 

d p {uj p+ )d p {uj p ~)d\ (w s+ )a s (a; s _)exp(-£Acj£)) 



-)- g 4M 2 



where Aw 



p+ 



(8) 



is the frequency 



detuning. L is the length of the medium, which comes 
from integrating z from to L. The integration of z does 
not involve the phases associated with the wavenumbcr 
of the modes because the sum of the wavenumbers is 
zero when dispersion is assumed to be negligible. Here 
X is again the interaction strength but incorporated 
with some constants from the electric field expressions. 
The frequency integrals have lower bounds extended 
from zero to negative infinity. This is mathematically 
legitimate because we are considering a system that 
operates at high frequency, where essentially there is no 
population present at low frequency. 

This interaction Hamiltonian docs not commute at 
different times, so strictly speaking, we should use the 
Dyson scries to calculate higher order terms in the 
unitary expansion. However, here we shall use the 
Taylor series as an approximation to the Dyson series. 
Later in this section, we shall show with the Taylor series 
that the non-instantaneous response of the Kerr medium 
will always produce spectral entanglement and poses a 
problem to the fidelity the gate. Because of this, for a 
non-instantaneous Kerr medium, we expect the more 
complicated Dyson series will also show a similar spectral 
entanglement problem. In the next section, we shall see 
that engineering the dispersion property of the medium 
can help eliminate the spectral entanglement and 



Effects of non-instantaneous response 

We have derived an expression for a Kerr medium with 
finite response time as shown in equation (|9|) . The finite 
response time translates into a spectral filtering function 
in frequency space. Here we shall see what happens when 
we apply this cross-Kerr process to two-singlc-photon in- 
put states. A single-photon state with a gaussian spectral 
profile is: 



ID = 



l 



gvItt J- 



dw^Hexp^-^ |0) (10) 



So the spectrally separable input state with modes p 
and s is: 



UD=- 



1 



<tV2it 



duip du! s al(u)p)al(u>s) 



— oo 
2 



2 \ 

-)|0) 



(11) 



where Vj = u> j — fi and [i is the mean frequency. Having 
the unitary acting on |11) gives: 

^ n )=E^b- l^)=y>n> (12) 



n=0 



n=0 



Note that |^o) = |H>- From equation © and ifH]). 
one can show that 



dujpdL0 s d p (uj p )al(L0 s ) 



CXP ^(2nM* + J) |0> (13) 
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Expression (|13|1 implies that cross-Kerr non-linearity 
inevitably induces spectral entanglement in the output 
state because the cross term in the argument of the 
exponential AoH^hP+o' 1 ) cannot vanish for non-zero 
A/. This is undesirable for quantum computation, 
as entanglement in the spectral domain will appear as 
decoherence in the logical basis, thus lowering the fidelity. 



Given the input state 
expected state is \ipex P ected) 
actual output state is: 



\ipout) = U\ip in ) 



= (|10) + |11))/V2, the 
(|10) - |11))/V2. The 




(14) 



Hence the expression for the fidelity is: 

2 



F J- 

4 



n\ \ ih J 

i-ixr v 



nM 2 + a 2 

2 



(15) 



In the case of 



F 



h 

1, equation (|15p becomes 
2 0. To see the effect 



where ?y = cr/M and X 
fast response, i.e. 77 <C 

V°° i (-iX) n -2- 

of the frequency response function on the state, we 
shall find the fidelity as a function of X and r\. In 
particular, we set r\ — 0.01 for Figure [21 Notice in 
the figure that even when X is large, the maximum 
fidelity is roughly in the order of 10 -4 . Similarly, for 
77 = 0.001, the maximum fidelity is on the order of 
10~ 6 for large X. Hence the fidelity is far too low 
for quantum computation in the fast response regime. 
The low fidelity implies that the output state is almost 
orthogonal to the expected state. It also implies that 
the fidelity of the output state with respect to the input 
state is near unity. Figure [3] plots 9 against X, where 
9 is the argument of the complex number (11|£/|11) in 
polar form 2l|. If the cross-Kerr medium does induce 
a 7r phase shift onto the two-singlc-photon state, then 
9 would be ±7r. However, the figure shows that for 
large X, 9 is roughly — -gfj, which indicates some but 
very little phase is induced onto the state by the medium. 
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FIG. 2: Fidelity vs X for 77 = 0.01. 
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FIG. 3: vs X for 77 = 0.01. Note that the argument 
is undetermined at X—0 because the complex number 

('ipexpectedllpout) = 0. 



that the higher order terms in the summation of equa- 
tion (|15[) is only significant up to the n c -th term. If we 
consider the regime where 77 ^> n c > 1, then we can see 
from equation (| 1 3|) that the cross term becomes insignif- 
icant, and F = \\l - YZ=o i~ iX T\ 2 ■ Setting X = tt 
gives F = 1, which means a tt phase shift is induced on 
1 11). Hence, we argue that it is the spectral entangle- 
ment caused by the non-instantaneous response function 
in the fast regime that leads to a low fidelity. Although 
this entanglement is reduced in the slow response regime, 
the presence of noise terms there reduce the fidelity of op- 
eration there as well, as pointed out by Shapiro [l2| . We 
seek a way to reduce the spectral entanglement whilst 
remaining in the fast response regime. 



Ill THE SPECTRAL RESPONSE OF A KERR 
MEDIUM WITH DISPERSION 



To understand the fact that it is the spectral entangle- 
ment that lowers the fidelity, we shall examine the case 
where the spectral entanglement is negligible. Which is 
the case of slow response, where M is small but con- 
veniently ignoring the noise for argument sake. Let us 
suppose that the Hamiltonian has the same form as equa- 
tion @, where the noise term is ignored, and also assume 



In the previous section, we have discussed what 
happens if we input a separable two-photon state into 
a Kerr medium with finite response time. In the fast 
response regime, the spectral entanglement induced by 
the non-instantaneous response function reduces the 
phase shift imparted and thus lowers the fidelity. Hence 
if we can engineer the phase such that the spectral 
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entanglement can be eliminated, then a unit fidelity may 
again be achieved. It is well known that by adjusting 
the dispersion in a parametric down conversion crystal, 
the s pec tral profile of the output photons can be mod- 
ified [17] . Here we propose to engineer the dispersion 
properties of a Kerr medium such that the resultant 
state is a separable two-photon state. Furthermore, with 
a special condition set for Xi we show that the resultant 
state can have a ir phase shift, which is what we hope for. 

The interaction Hamiltonian with dispersion and non- 
instantaneous response is: 



H(t) ~xL J J J J duj p +dujp-d(jj s +duj s - 
at p {uj p+ ) a p (ujp- )at s (uj s+ ) a s ) 

p[ up ) + cp[ UP 

( LAk\ ( iLAk\ , . . , . 

exp [ — - — ] exp(— iAut) (16) 



smc 



where Ak ~ k p {ui p + — + k s (u) s + ~ w s _) is the differ- 
ence in the wavenumber of mode p and s, i.e. the phase 
mismatch, up to the first order terms after a Taylor se- 
ries expansion about the mean frequency. Here k'j is the 
derivative of the wavenumber of mode j with respect to 
u)j and evaluated at the mean frequency. We have as- 
sumed that the higher order terms of Ak are negligible. 
Moreover, the zero-th order term is zero because of con- 
servation of momentum. Later in the appendix, we shall 
show that when the spectral effects of the dispersion dom- 
inates over the non- instantaneous effect, the interaction 
Hamiltonian does commute at different times, and thus 
we may use the Taylor series to calculate higher order 
terms of the unitary expansion. Hence the n-th state in 
the sum now has the expression: 



that in the limit of £ — > oo, the cross term C3 is zero 
and thus the resultant state is separable. Along with a 
special condition on \ that allows the unitary to give a 7r 
phase shift, we could achieve unit fidelity. It is therefore 
of interest to know what length of material can give a 
high fidelity. The new expression for the fidelity is: 



51 ~^~\ 



2-kixL 



exp 



fM 2 "7"o- 2 (l + L 2 M 2 1 {k p - k's) 2 ) 1 " 
APn + cr 2 + L 2 M 2 7cr 2 (fc; - fc;) 2 

nL 2 M 2 a 2 {k' p - k' s ) 2 



4(M 2 n + a 2 + L 2 M 2 ja 2 {k p - k' s ) 2 ) 



We now set the special condition 



2X 



(18) 

Wp-K\ exp (*r)' 
which is derived from taking the limit L — > 00, such 
that the unitary induces a ir phase shift on |11), in other 
words, F = 1 in equation (fl8|) . With this condition, we 
simplify the expression for the fidelity to: 



^ (-nr)«C 

> ■ exp 



47(n + £ e ) 



(19) 



where we have assumed <C L 2j y(k p — k s ) 2 , i.e. the 
effect of dispersion dominates over the effect of non- 
instantaneous response. Figure [4] is the plot of fidelity 
versus the dimensionless C. When F = 0.999, we have 



7< 



C = 100, so for the reasonable parameters 

10 -8 (s/m) and a = 10 13 (iJz), the corresponding length 
of the medium is L = 2.28 * 10~ 3 (to). 



IVO =■ 



M 2n £ 2n (a 2 + M 2 C 2 ) 



I {a 2 + M 2 C 2 ) n (2nM 

1 /7T\" 1 



M 2 C 2 ) 



dwpdojs at ( uj p ) d\ ( us s ) 



{C\Vp+Clvl + C2V p +C2Va+C3VpV3 + Ci) |Q\ (17) 



where we define C = Lu^J j(k' p — k' s ) 2 , 

r _ nM 2 +cr 2 +M 2 C 2 (~l _ inCM 

01 ^ 4a 2 (2nM 2 +a 2 + M 2 C 2 ) ' ° 2 



2cr x /7(2nJ\/ 2 +(T 2 +A/ 2 £ 2 ) 
n(2nM 2 +cr 2 ) 



/~t _ nM 2 j (~i _ _ 

L '3 — 2<r 2 (2nM 2 +<j 2 +A/ 2 £ 2 ) ' aUU 04 — 4 7 (2nA/ 2 +a 2 +A/ 2 £ 2 ) 

Here we have approximated the sine function of the 
phase matching function as a gaussian function via the 
approximation sinc{x) ~ exp(— 7a; 2 ) with 7 « 0.193... 
calculated from matching the full-width-half-maximum 
of the two functions. From equation (|17[) . we can see 




1000 100000. 



FIG. 4: Fidelity vs £. 



IV. CONCLUSION 



In this paper, we have modeled the spectral effects of 
a non-instantaneous cross-Kcrr medium and discussed 
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how the noise terms of the interaction Hamiltonian 
can be ignored in the fast response regime. We have 
shown that for a separable input of two-single-photon 
state, a fast response cross-Kerr medium inevitably 
spectrally entangles the state at output. Ideally, a 
cross-Kerr medium would impose a 7r phase shift on 
the input state, but we have shown that the presence 
of the entanglement reduces the phase shift to very 
little, and found that the resulting fidelity is too low 
for quantum computation. We have illustrated that 
a way to circumvent the problem is to engineer the 
dispersion properties of a fast response medium to give 
a dominant spectral effect, in addition to satisfying 

the phase matching condition La^J^f(k' p — k' s ) 2 — ► oo. 
The dispersion helps the state to remain separable, and 
by further picking a particular value of Xj a n phase 
shift can be induced on the two-single-photon state. 
We therefore argue that it is possible to build a high 
fidelity cross-Kerr quantum gate from a single pass 
medium of sufficient non-linearity, provided one could 
engineer the dispersion of a fast response medium. Our 
calculation shows that in order to achieve a fidelity of 
0.999, we need the dimensionless parameter C = 100. 
An example is where the difference of the reciprocal of 
the group velocities is 10~ 8 (s/m) and the spectral width 
of the photons is 10 13 {Hz), and a length of medium of 
2.28(mra). 
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APPENDIX 

Al: Approximating the Dyson Series with Taylor 
Series for the Interaction Hamiltonian with 
Non-Instantaneous Response 

The interaction Hamiltonian with a non-instantaneous 
response does not commute at different times and 
strictly speaking, we should use the Dyson series instead 
of the Taylor series to calculate higher order terms of 
the unitary expansion. However, we have calculated the 
fidelity for second order term of the Dyson and Taylor 
series with various values of M and a of interest and the 
fidelity seems to be close to unity. In particular, when 



M = 10 10 and a = 10 9 , the fidelity is approximately 
1.046 with a roughly estimated numerical error of 
±0.077. Fidelity is always bounded by one and this 
bound is within the numerical error estimated. Hence 
arguably, the second order term of the Taylor and Dyson 
series are identical. This suggests that the Taylor series 
may possibly give a good approximation for higher order 
terms of the evolution of spectrally-Gaussian input 
states passing through a bulk cross Kerr medium. 

Furthermore, if the Kerr medium consists of thin slices 
separated sufficiently apart, such that the wavepacket ex- 
its one slice before entering another, and that each slice 
only provides a weak interaction, then it is appropriate 
to apply the Taylor series instead of the Dyson series 
for the unitary operation expansion. This is because the 
unitary operator can be expressed as equation (f2"0"|) and 
if the slices are sufficiently apart, the time bounds of the 
integrals can be conveniently extended to — oo and oo 
by considering far field limits [HI, and thus each factor 
in the series is the same, and the time ordering opera- 
tor can again be dropped, which gives the Taylor series 
expansion. 

U(h,t ) = lim r{ ( 1 - % - f 1 H(t)dt) 

n— oo I \ h J tl _ AT J 

r AT H(t)dt)... 

\ n Jt 1 -2AT J 

(i-i£ +AT ^} (20) 

where AT = (ii - t )/n. 



A2: Time Commutivity of the Interaction 
Hamiltonian with Dispersion Effects 

Here we shall prove that in the fast response limit with 
a dominant effect from the dispersion of the medium, the 
Hamiltonian in equation p6[) commutes at different times 
and the expression for the unitary operator is properly 
given in equation @ without the time ordering operator. 

In the limit of very large M, the interaction Hamilto- 
nian in equation (|16[) tends to: 

H{t) JJJJduj p+ dujp^dbj s+ dLd s -dl(LUp + )a p (LLjp-) 

at(^ + )a s ( Ws _)sinc(^)e^e- lA - t (21) 

Hence the commutation of the interaction Hamiltonian 
at times t\ and ti is: 
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f f 

=X 2 L 2 / ... / dtJp+dtJp-dLLj s+ did s _duj' p+ dw' p _du)' s+ dLu' s _ 
al {u!p+)d p (ujp_)d\{LJ s+ ) d s (w a _ ) , 
al(uJp + )a p (w p _ )a\ (uj' s+ )a s (w' s _ ) 



u' p _), Ak' s , = k'M+ -<-)■ Now 



where Ak' p = 
Ak' p , = k' p (uj' p+ 
integrating over lo' s+ and w s _ for the first two terms gives 
zero, as they cancel out each other as well. Similarly by 
changing variables and integrating other terms will have 
the terms cancel out each other. Hence the interaction 
Hamiltonian with a dominating spectral effect from the 
dispersion of the medium, commutes at different times. 



/ LAk\ [LAk'\ :LAk f LAk' 
sine ^ j sine ^ — - — Je 2 e 2 

where Ah' w k' p {u' p+ - u' p _) + k' s (cu' s+ - lu' s _) and 
Au' = uj' p+ - J p _ + uj' s+ - u' s _. 



(22) 



The commutation relation in the above equation can 
be expressed as a sum of 8 terms: 

d p {oj p+ )a p (u p _ )a\ (w s+ )a s (w s _), 
a\ {u'p + )d p {uj' p _) d\ (w' s+ ) d s (u' s _ ) 

= - ( d %+. s '-al( w U)^(^-)«p( w P+)apH-) a pH+) a pH- 

- 5 s / +!S _a t s (^ s+ )a s (w^)aj ) (w p+ )ap(wp_)aj J (^ + )a p (^_) 
s p+,p'- «p H+)«p )°>l (us+jds (w s _ )a* (w 8+ )a a {u>' 8 _ ) 

- H- ) & t ( UJ s+)d s (uj s -)al (u)' s+ )a s (J s _ ) 
6 p'+,P-ai(u p+ )a p (u; p ^)5 s , + ^dl(uj s+ )d s (u;' s _) 

-S p , +jP _al(uj p+ )a p (uj' p _)S s+!S ,_al(uj' s+ )a s (uj s ^ 
^+. P '- a j,K+)a P K-)^+, s '-fflJ( w l+)fl s ( w s -) 

-*p+,p'-OpH+) a p( a; p-)^'+ 1 *-aU w «+)o*( w *-)) 

(23) 

where 8 X%V = 5{uj x — uj y ). After substituting this into 
equation (|22|) . we will also obtain 8 terms. If we inter- 
change the variables u>' s , «-> cj s + and u)' s _ <-» w s _ for the 
first term and then add with the second term, we obtain: 



[H(h),H(t 2 )] 



du) p+ duj p -duj s+ doj s ^duj' p+ duj' p _dijj' s+ dijj' s _ 



a s (w s+ )a s (^_)aJ(w p+ )a J) (wp_)a|,(^ + )ap(^_)5 s / + , s _ 

sinc( L(A ^ 2 +Afc - ) )e^^^) 
sinc( L(Afc ^ +M -' ) )e^(^^) 

_ -i(w p+ -Wp_+w s + -w s _)ti ,— a)'_ +w! j, -w^_)t 2 
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